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ABSTRACT 


The  influences  of  release  of  latent  heat  of  vaporization  on 
vertical  motion,  production  of  kinetic  energy,  and  development  of  low-level 
vorticity  have  been  investigated .  The  vertical  motion  was  obtained  by  solving 
the  customary  (diagnostic)  omega -equat  ion,  making  allowance  for  horizontal 
variations  of  the  static  stability.  The  results  were  compared  with  vertical 
velocities  computed  by  kinematical  techniques,  and  it  was  found  that  it  was 
necessary  to  include  the  effect  of  released  latent  heat  in  order  to  obtain 
satisfactory  agreement  between  the  two  sets  of  data.  Experimental  com¬ 
putations  showed  that  the  effects  of  spatial  variations  of  the  static  stability 
were  important  outside  regions  of  condensation. 

Numerical  solutions,  with  and  without  inclusion  of  released 
latent  heat,  were  used  to  obtain  ageostrophic  wind  components  and  their  effect 
upon  the  production  of  kinetic  energy  and  vorticity.  The  following  results 
emerged:  (a)  the  influence  of  released  latent  heat  was  of  the  same  order  of 
magnitude  as  the  effect  of  dry-adiabatic  circulations;  (b)  to  the  amplification 
of  the  vertical  motions  that  resulted  from  released  latent  heat  corresponded 
intensification  of  the  low-level  convergence  and  high-level  divergence;  (c) 
the  ageostrophic  winds  associated  with  these  fields  of  convergence  and 
divergence  had  components  toward  lower  pressure,  thus  giving  positive 
contributions  to  the  production  of  kinetic  energy  both  at  low  and  at  high  levels; 
(d)  the  computed  rate  of  production  of  low-level  vorticity  exceeded  the  observed 
rate,  and  evidence  suggests  that  frictional  effects  may  be  important. 
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1. 


INTRODUCTION 


During  the  Nineteenth  Century  many  investigators  (  e.  g.  Reye 
1872,  and  Ferrel  1889)  maintained  that  the  release  of  latent  heat  of  vaporization 
contributed  substantially  to  the  development  of  cyclones,  and  some  support  for 
this  thought  was  provided  later  when  it  became  evident  that  the  rate  of  release 
of  latent  heat  far  exceeded  the  production  of  kinetic  energy ,  On  the  other  hand, 
Margules  (1903)  calculated  that,  in  a  closed  system,  the  increase  in  kinetic 
energy,  resulting  from  the  redistribution  of  juxtaposed  warm  and  cold  air 
masses,  was  negligibly  affected  by  the  occurrence  of  condensation  in  the  warm 
air,  provided  the  latter  were  potentially  stable.  It  should  be  mentioned,  how¬ 
ever,  that  Margules  considered  only  initial  and  final  states  and  made  no 
inference  as  to  the  rate  of  increase  in  kinetic  energy.  Furthermore,  the 
assumption  of  closed  boundaries  would  limit  the  application  of  his  results 
to  natural  systems .  It  is,  perhaps,  due  to  Margules'  influence  that,  in  later 
years,  little  attention  was  paid  to  the  importance  of  release  of  latent  heat  as 
source  of  energy  in  cyclone  development .  Only  recently,  mainly  as  a  result 
of  advances  in  numerical  forecasting,  has  interest  again  become  centered  on 
the  role  of  latent  heat  in  the  development  of  motion  systems . 

The  purpose  of  this  paper  is  to  explore  the  manner  in  which 
release  of  latent  heat  affected  the  development  of  a  selected  cyclone .  More 
specifically,  we  shall  be  concerned  with  the  effects  on  the  fields  on  vertical 
motion  and  divergence  and  the  rates  of  production  of  kinetic  energy  and  low- 
level  vorticity. 


-1- 


2. 


SOME  EARLIER  INVESTIGATIONS 


An  early  idea  which  has  recurred  frequently  in  the  literature  is 
that  the  release  of  latent  heat  of  vaporization  tends  to  decrease  the  density 
of  the  air  column  in  which  condensation  occurs  and  thereby  lower  the 

pressure  at  the  base  of  the  column .  Reye  (1865)  considered  an  atmospheric 

2  o 

column  with  cross  section  of  1  m  in  which  1  kg  of  water  vapor  at  10  C 

was  condensed,  thereby  producing  1  mm  of  precipitation.  From  the  released 
latent  heat  he  computed  the  horizontal  expansion  and  found  a  pressure  reduction 
of  about  0.8  mb.  This  provided  him  with  a  relation  between  precipitation 
and  surface  pressure  from  which  he  concluded  that  the  heat  released  by  con¬ 
densation  was  the  chief  cause  of  cyclone  development.  In  a  similar  manner, 
Refsdal  (1930)  considered  a  saturated  atmosphere  with  unstable  stratification 
above  the  800  mb  level .  By  letting  a  parcel  of  air  rise  from  800  to  500  mb 
Refsdal  found  that  a  saturated -adiabatic  lapse  rate  would  become  established 
in  the  ascending  column  and  that,  under  such  temperature  conditions  as 
normally  prevail  in  middle  latitudes,  about  10  mm  of  precipitation  would  have 
been  released .  Refsdal  assumed  that  the  rising  air  would  be  about  5°  C  warmer 
than  its  environment  between  750  and  500  mb,  and  his  calculations  showed  that 
if  the  height  of  the  500  mb  surface  remained  constant  the  sea  level  pressure 
would  be  reduced  by  about  4.6  mb.  like  Reye,  Refsdal  considered  that  he  had 
obtained  a  valid  relationship  between  the  intensity  of  precipitation  and  the 
rate  of  fall  of  the  sea-level  pressure. 

On  the  other  hand,  Hann  (1874)  studied  records  from  Batavia  and 
found  that  heavy  rains  (  20  mm  or  more  per  hour)  were  accompanied  by 
rising  pressure,  and  from  this  he  concluded  that  the  release  of  latent  heat 
could  not  account  for  the  development  of  low-pressure  systems.  It  may  be 
noted,  however,  that  the  data  examined  by  Hann  likely  reflected  the  effects 
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of  downdrafts  from  convective  clouds . 

Fundamental  progress  in  the  understanding  of  the  energy 
sources  of  cyclones  was  achieved  by  Margules  (1903),  whose  influence  is 
clearly  evident  in  the  later  works  of  the  Norwegian  school  of  meteorologists . 
Margules  considered  first  a  double  chamber  (as  dhown  on  the  left  in 
Fig.  2.2.  1)  in  which  two  masses  of  dry  air  in  stable  or  neutral  equilibrium, 
having  the  temperatures  and  T0  (with  <  T0  ), 

were  separated  by  a  vertical  wall .  It  was  assumed  that  the  pressure  remained 
constant  along  the  top  surface  and  that  the  discontinuity  in  pressure  at  the 
partition  increased  from  zero  at  the  top  to  a  maximum  value  at  the  bottom  of 
the  chamber .  If  the  effects  of  friction  are  neglected  and  the  boundaries  assumed 
to  be  rigid,  it  is  readily  shown  that 

(2.2.1)  K  +  P+  T-  Q*  constant 

where  K  »  P  ,  and  X  denote  the  kinetic,  potential,  and  internal  energy, 
respectively,  and  Q  is  the  heat  added  to  the  system  . 
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Fig.  2.2.1.  Margules' two-chamber  model. 
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When  the  partition  between  the  chambers  is  removed  the  air  masses 
will  become  redistributed  as  shown  in  the  diagram  on  the  right.  If  subscript 
(X  identifies  an  adiabatic  process,  the  difference  between  the  initial  and 
the  final  states  is  found  to  be 

(2.2.2)  AKa  -  ^  + 


Furthermore,  if  the  two  chambers  are  of  equal  size,  the  following  approximate 
relationship  can  be  shown  to  hold 


(2.2.3) 


tb-ta 

Ta 


where  M  is  the  total  mass. 

Margules  next  considered  the  effect  when  the  mass  B  was  re¬ 
placed  by  saturated  air  having  a  saturated -adiabatic  lapse  rate.  When  the 
partition  is  now  removed,  condensation  will  take  place  in  the  mass  0  , 
and  Eq.  (2.2.1)  gives 


(2.2.4) 


AKh  =  -A(p  +  DK  +  a 


Here,  subscript  h  signifies  that  the  release  of  latent  heat  is  involved. 
Margules  found  that 


(2.2.5) 


-j>(P  +  I)a- A  (P+I)h]  =  Q 
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showing  that  when  condensation  occurs,  the  decrease  in  potential  and 
internal  energy  is  less  than  it  would  have  been  had  the  process  been  dry- 
adiabatic,  the  difference  being  equal  to  the  amount  of  released  latent  heat. 
Equation  (2.2.5)  may  be  interpreted  to  state  that  the  release  of  latent 
heat  of  vaporization  merely  diminishes  the  rate  of  cooling  of  the  ascending 
mass  and  does  not  contribute  to  the  increase  of  kinetic  energy.  Thus 

(2.2.6)  AKh=  A  K  a 

a  result  which  might  seem  to  defeat  the  main  purpose  of  the  present  paper. 

It  should  be  noted,  however,  that  Margules  did  not  consider  the  influence  of 
the  surroundings  on  the  system,  nor  did  he  draw  any  conclusions  about  the 
rate  at  which  the  conversion  of  energy  takes  place;  it  is  possible,  and  indeed 
probable,  that  the  amplifying  effect  that  release  of  latent  heat  has  on  the 
vertical  motion  may  accelerate  the  conversion  process.  We  shall  return 
to  the  discussion  of  this  in  later  sections. 

In  a  second  paper,  Margules  (1906)  considered  the  case  when 
the  mass  B  is  potentially  unstable,  and  found  a  significant  contribution 
from  the  released  latent  heat.  In  the  following  Sections  3  and  4,  discussion 
will  be  restricted  to  cases  where  the  air  is  potentially  stable,  and  it  will  be 
shown  that  even  in  such  cases  release  of  latent  heat  may  lead  to  significant 
increases  in  the  kinetic  energy  of  mid -latitude  cyclones. 

In  connection  with  numerical  predictions  the  vertical  velocity  is 
usually  obtained  by  solving  the  so-called  -equation  derived  by  Fjbrtoft 
(1955);  a  form  of  this  equation  will  be  reproduced  in  Section  3.1 .  In  the  early 
applications  of  Fjbrtoft' s  equation  it  was  customary  to  ignore  the  effects  of 
released  latent  heat.  It  soon  became  evident,  however,  that  the  neglect  of 
latent  heat  in  the  UJ  -equation  resulted  in  a  substantial  underestimate  of  the 
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vertical  motions  within  the  regions  where  condensation  occurred .  The  first 
attempt  at  accounting  for  the  effect  of  condensation  was  made  by  Smagorinsky 
(1956)  who  found  that  the  calculated  velocities  could  be  augmented  by  an  order 
of  magnitude,  partly  by  accounting  for  the  horizontal  variation  in  the  static 
stability  but  primarily  by  including  the  effect  of  release  of  latent  heat .  Other 
authors  have  also  demonstrated  the  importance  of  released  heat  of  vaporization 
in  dynamical  processes  .  For  example,  Manabe  (1956)  studied  a  case  of  heavy 
large-scale  precipitation  in  which  a  considerable  increase  of  potential  vorticity 
occurred  in  the  lower  troposphere  .  He  found  that  the  increase  in  potential 
vorticity  predicted  from  heat  released  by  condensation  was  of  the  right  order 
of  magnitude  to  explain  the  observed  changes  .  Furthermore,  Aubert  (1957) 
found  that  the  released  heat  of  vaporization  greatly  increased  the  upward  motion, 
and  he  found  that  this  release  also  tended  to  lower  the  heights  of  the  isobaric 
surfaces  in  the  lower  troposphere  and  to  raise  them  in  the  upper  atmosphere . 
These  changes  resulted  in  deepening  of  low-level  cyclones  and  acceleration  of 
their  rates  of  movement.  Smebye  (1958)  conducted  a  number  of  experiments  in 
quantitative  precipitation  forecasting;  he  included  the  effect  of  heat  released 
by  condensation  and,  except  in  the  vicinity  of  convective  systems,  found  good 
agreement  between  predicted  and  observed  amounts  of  precipitation.  Further 
refinements  were  introduced  by  Pedersen  (1962)  who,  like  Smebye,  demonstrated 
the  importance  of  the  inclusion  of  released  latent  heat. 

In  a  recent  study  of  the  effects  of  heat  sources  on  the  behavior  of 
winter  cyclones  on  the  North  Atlantic  Ocean,  Petterssen  and  collaborators 
(1962)  computed  the  rate  of  change  of  the  thickness  from  1000  to  500  mb  during 
typical  stages  of  development  of  extra -tropical  cyclones,  and  found  that  the 
inclusion  of  released  latent  heat  as  well  as  the  eddy  flux  of  sensible  heat  were 
important.  In  typical  cases,  the  maximum  flux  of  sensible  heat  occurred  in 
the  rear  of  the  cyclone  center  and  amounted  to  about  1  cal  cm  min'  .  Since 
a  rate  of  precipitation  of  1  mm  hr  1  would  liberate  an  equivalent  amount  of  heat, 
it  is  evident  that,  when  the  temperature  is  high,  the  release  of  latent  heat  may  be 
of  paramount  importance . 


3. 


THEORY  AND  PROCEDURE 


The  purpose  of  this  section  is  to  present  the  basic  theory  and  an 
outline  of  the  procedure  to  be  employed .  Further  details  of  the  methods  used 
will  be  discussed  in  Section  4 . 

3 . 1  Vertical  Velocity 

With  the  symbols  defined  elsewhere,  the  quasi -geostrophic 
vorticity  equation  and  the  corresponding  thermodynamic  energy  equation  may 
be  written: 


(3.1.1) 


and 


(3.1.2) 


R_H 

CpP 


It  may  be  noted  that  Eq  .  (3.1.1)  satisfies  the  consistency  re¬ 
quirements  of  Hollmann  (1956),  Wiin-Nielsen  (1959),  and  Lorenz  (1960^,  Uctte/ 
provided  variations  in  Cj^  and  arising  from  changes  in  + 

However,  Eq.  (3.1.2)  does  not  satisfy  these  criteria,  unless  for  each 
isobaric  surface  6  is  taken  to  be  a  constant,  and  H  is  defined  such  that 


JhWA  =  o 


The  above  integral  is  to  be  taken  over  the  area  in  question.  Nevertheless, 
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as  will  be  shown  in  Section  4,  the  use  of  an  inconsistent  form  of  Eq.  (3.1.2)  re¬ 
presents  the  instantaneous  physical  relations  better  than  does  the  consistent  form, 
although  the  use  of  the  former  in  numerical  forecasting  may  cause  spurious  energy 
changes  after  a  large  number  of  time  steps . 

Combining  Eqs.  (3.3.1)  and  (3.1.2),  we  obtain 

(3 . 1 .3)  4- V1 JB- ' 

The  above  equation  represents  a  simplification  of  the  Ui  -equation  derived  by 
Fjbrtoft  (1955). 

Since  Eq.  (3.1.3)  is  linear  in  Ld  ,  we  may  write 

U>  -  U>tt  + 

where  U)a  may  be  interpreted  as  the  vertical  velocity  due  to  dry-adiabatic 
processes,  and  as  the  additional  velocity  caused  by  non-adiabatic 

heating.  Equation  (3 . 1 .3)  may  then  be  replaced  by  the  following  two  equations: 
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In  Sections  4 .3  and  4.4,  the  above  two  equations  will  be  solved  by 
relaxation  techniques,  using  a  vertical  mesh  size  of  200  mb  and  keeping  the 
dependent  variable  zero  at  p  »  O  and  at  jps  1000  mb  •  The 
two  solutions  will  then  be  added  to  give  U)  .  For  comparison,  the  vertical 
velocity  will  also  be  determined  kinematically,  taking  orographic  influences 
into  account.  Ageostrophic,  frictional  and  orographic  effects  neglected  in 
the  numerical  solution  and  errors  inherent  in  the  kinematic  technique  will  be 
the  main  factors  preventing  an  exact  agreement  between  the  two  results. 

The  main  contribution  to  H  will  be  assumed  to  be  from  heat 
released  by  condensation,  and  other  heat  sources  and  sinks  will  be  neglected. 
Furthermore,  it  will  be  assumed  that  condensation  occurs  as  liquid  water  which 
immediately  falls  out  as  precipitation.  Apart  from  the  question  of  the  phase  of 
the  condensed  water  substance,  this  last  assumption  is  at  variance  with  the 
conclusions  of  Cunningham  (1952).  He  found  that  in  cases  of  steady  rain, 
appreciable  amounts  of  water  are  collected  when  precipitation  elements  pro¬ 
duced  by  clouds  aloft  fall  through  low-level  stratus  and  scud.  These  low  clouds 
owe  their  existence  not  to  large  -scale  vertical  motion  but  rather  to  evaporation 
from  falling  rain  and  from  the  earth's  surface.  However,  if  the  former  moisture 
source  is  more  important  than  the  latter,  Cunningham's  results  should  not 
seriously  affect  the  validity  of  our  assumption.  This  is  because  the  loss  of 
water  by  evaporation  of  falling  rain  (  which  has  been  neglected)  will  be  com¬ 
pensated  by  a  gain  when  the  raindrops  fall  through  the  low  clouds. 

To  determine  the  spatial  distribution  of  H  ,  use  will  be  made  of  the 
vertical  velocity  as  computed  kinematically  and  the  observed  precipitation 
rates,  both  averaged  over  areas  centered  at  grid  points  (see  Fig.  4.3.1). 

It  will  be  assumed  that  condensation  occurs  only  in  the  region  between  900  and 
300  mb.  The  atmosphere  will  be  divided  into  layers  200  mb  dpep  centered  at 
800,  600  and  400  mb,  and  the  moisture  data  from  radio  soundings  will  be  used 
to  determine  which  of  these  layers  could  readily  produce  precipitation.  For 
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convenience,  a  possible  precipitation-producing  layer  will  be  defined  as  one  which 
contains  sub-layers  where  T  -  T  *  i  S°C  .  This  figure  has  been 
chosen  somewhat  arbitrarily  to  tnake  allowance  for  the  fact  that  100  per  cent 
relative  humidity  is  rarely  reported  by  radiosondes  even  in  cases  of  continuous 

i  I 

rain.  Let  the  subscripts  i  and  m  denote  conditions  at  the  i/  grid 
point  and  m  ^  level  (  m  =  1,  2,  3  correspond  to  800,  600  and  400  mb 
respectively) .  At  each  grid  point,  we  have 

3 

(3-1.6)  -  £l  I  Cim(T,p)<A).m^=  Pi 

■m  -  \ 

where  Pj  is  the  total  precipitation  and  is  the 

precipitation  produced  by  a  200  mb  deep  saturated  layer  ascending  at  the  rate 
-3  - 1 

of  10  mb  sec  (obtained  by  modifying  values  given  by  Fulks,  1935).  The 
quantity  is  defined  alternatively  as  follows: 

<  -th 

Oy„  =  1  if  the  vrv  layer  has  sufficient  moisture  to  produce 
precipitation  and  <  O  ;  otherwise 

ivw  =  O 

Since  LJ,;m  is  known  from  kinematic  computations,  Eq .  (  3 . 1 .6)  can  be 
solved  for  £  L  ,  which  is  the  fraction  cf  the  column  that  is  producing  pre¬ 
cipitation  (  if  i  -  £,  =  i,  £t=l  indicates 

saturation  from  900  to  300  mb) .  Having  found  £^  ,  is  determined 

from 

Him  =  -  — -  C;„(T,p)UimLv  C 

$0£>  tnb 


(3.1.7) 
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The  quantity  200  tnb  /  ^  is  the  mass  of  a  column  of  air  of  unit  area  and  200 
mb  depth  so  is  the  rate  of  heating  per  unit  mass .  Hence,  Eqs.  (3.1.3) 

or  (3 . 1 .5)  can  now  be  solved . 

An  alternative  way  to  solve  Eq  .  (3.1.3)  is  the  following.  If 
l*>Cn\  in  Eq.  (  3.1.7)  is  interpreted  as  the  solution  of  Eq.  (  3.1.3)  rather 
than  the  kinematically  determined  vertical  velocity,  Eq .  (3.1 .3)  may  be  written 

(3.1.8)  +  + 

where 

(3.1.9)  <fs  s  A-  A -  £vCirn( T,p)Lv^ 

Cp  “200  *nV> 


There  are  two  difficulties  in  the  way  of  obtaining  solutions  to  Eq . 
(3.1.8)  rather  than  Eqs  .  (  3 . 1 .3)  or  (3 . 1 .5) .  In  the  first  place,  the  unknown 
cannot  be  determined  until  Eq  .  (3.1.8)  has  been  solved  for  U>  .  If 
we  make  an  initial  guess  of  £  ^  ,  Eq  .  (  3 . 1 .6)  will  not  be  satisfied  unless 

our  initial  guess  were  correct .  Secondly,  Eq  .  (  3.1.8)  will  not  have  a  bounded 
solution  if  <S  -  ffs  ^  o  .  Furthermore,  since  numerical  methods 
will  be  used,  large  errors  may  result  if  O  -  Gs  ,  though  positive,  is 
close  to  zero.  Of  course,  if  a  -  <JS  <  0  ,  the  atmosphere  is  unstable  and 
the  pattern  of  vertical  motion  takes  the  form  of  small-scale  convective  cells. 

In  this  case,  the  solutions  of  Eqs .  (3.1.3),  (3 . 1 .4)  or  (3 . 1 .5)  may  have 
little  meaning.  However,  instability  should  manifest  itself  in  the  precipitation 
pattern,  and  if  this  pattern  is  fairly  uniform,  it  may  be  concluded  that  thej  solutions 
of  Eqs.  (3.1.3),  (3.1.4)  or  (3.1.5)  then  represent  real  atmospheric  processes. 
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3.2  Changes  in  Kinetic  Energy 

Following  Phillips  (1956),  we  multiply  Eq.  (3.1.1)  by  -■jZ/f  . 
neglect  small  terms  arising  from  the  variations  of  •f  ,  and  obtain  the  quasi - 
geostrophic  energy  equation 

(3.2.1)  r  -V.fe4V<  -  j-  -  CO* 

it  3  3 

where  *  (uy  +  )/2  is  the  kinetic  energy  per  unit  mass  of  the 

geostrophic  motion. 

On  the  other  hand,  multiplying  scalarly  by  'V  the  equation  of 
frictionless  motion 


(3.2.2)  dV  =  -V$  -fJtxV 

dt 

one  obtains 

(3.2.3)  ^  =-V'(^+  $)V-|- 
it 

where  fe  =  (u\v1)/ 2  is  the  kinetic  energy  of  the  horizontal  motion 

per  unit  mass . 


The  wind  field  will  be  assumed  to  be  composed  of  a  geostrophic 

component  and  a  small  irrotational  perturbation  *  given  by 

$ 


(3.2.4) 


v'=  -  vcp 
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where  the  velocity  potential  (p  satisfies  the  condition 

-.t,*  ^  <>3 

0.2.5)  rj  =  — 


It  has  been  postulated  here  that  the  geostrophic  wind  contains  all  the  vorticity, 
and  Petterssen  (1953)  has  shown  that,  on  occasions,  this  is  not  justified.  How¬ 
ever,  this  assumption  has  been  made  in  order  to  obtain  a  convenient  separation 
of  the  motion  into  geostrophic  and  ageostrophic  components.  Moreover,  the 
effect  discussed  by  Petterssen  is  important  only  in  cases  where  the  deformation 
is  strong  and  the  vorticity  weak . 

Substituting  'V  =  into  Eq .  (3 . 2 . 2)  and 

proceeding  as  before,  we  obtain  the  approximate  energy  equation 

(3.2.6)  ^ 

d±  3  ’  3 

which,  by  virtue  of  its  similarity  to  Eq.  (3.2.3),  will  represent  the  physical 
conditions  better  than  does  Eq  .  (3.2.1). 

An  alternate  expression  to  Eq.  (3.2.6)  is 

(3.2.7)  -V-V# 

it  *  5  ’ 

If  Eq.  (3.2.7)  is  integrated  over  a  region  bounded  at  the  top  and  bottom  by 
the  pressure  levels  b,  and  ,  respectively,  we  obtain 
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(3.2.8) 


SK« 

st 


where 


f  f  v'-vf  sa 


i 


is  the  total  kinetic  energy  of  the  geostrophic  motion  within  the  region,  and  \^n 
is  the  component  of  ^  along  the  exterior  normal  to  the  lateral  boundary  L.  • 


Equation  (3.2.8)  thus  gives  the  rate  of  change  of  the  kinetic  energy  of 
the  geostrophic  motion  within  the  region .  The  first  two  integrals  on  the  right 
represent  the  net  horizontal  and  vertical  import  of  kinetic  energy  into  the  region. 
The  last  term  is  the  rate  of  working  by  the  pressure  forces  within  the  region . 


Equation  (3.2.8)  may  now  be  divided  into  three  parts:  one  which  is 
due  to  purely  geostrophic  motion;  one  which  is  due  to  ageostrophic  effects  arising 
from  dry-adiabatic  processes;  and  one  which  is  caused  by  released  latent  heat. 
From  the  solutions  of  Eqs  .(3.1 .4)  and  (3.1.5)  and  Eqs .  (3 .2 .4)  and  (3 .2 .5), 
we  find  and  ,  the  ageostrophic  wind  components  caused  by  dry - 

adiabatic  processes  and  released  heat  of  vaporization,  respectively.  The  three 
components  of  Eq.  (3.2.8),  denoted  by  the  subscripts  ^  ,  Cl  and  h  , 

respectively,  may  then  be  written: 


(3.2.9) 
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(3.2.10) 


(3.2.11) 


_jjv;.v§sAS> 


An  equivalent  though  less  simple  approach  would  be  to  integrate 


Eq .  (3.2.6).  Analogous  to  Eq .  (3.2.8)  we  would  then  obtain 

^  t 

(3.2.12, 


where 


VJ 


is  the  normal  component  of 


The  third  and  fourth  integrals  on  the  right  side  of  Eq .  (3.2.12) 
represent  the  rate  of  working  by  the  pressure  forces  on  the  boundaries  of  the 
region.  It  is  readily  shown  (  e.g.,  White  and  Saltzman,  1956)  that  the  last 
term  is  the  rate  of  conversion  of  potential  and  internal  energy  into  kinetic 
energy . 


Expressions  equivalent  to  Eqs.  (3.2.10)  and  (3.2.11)  are 
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[«•]  .  [[<^-<^1  j  -  f  R  *  SL  7 

4<*hV  7  -  f  R SA  J 

[^]  --R^V^Uj -f#y^ 

®  K  - 

(3.2.14) 

+|[<i-JPt-(^VJ7 

The  chief  difficulty  in  the  way  of  using  Eqs.  (3.2.13)  and  (3.2.14) 
rather  than  Eqs.  (3.2.10)  and  (3 .2, 11),  is  encountered  in  calculating  the  integrals 

§41  If 

-JR*a7 

Since  these  terms  depend  only  upon  the  values  of  ^  and  on  the 

borders  of  the  region,  it  is  likely  that  the  boundary  error  in  the  numerical 
solution  of  these  Induced  velocities  will  adversely  affect  the  results.  Con¬ 
sequently,  in  Section  4 .5  the  terms  in  the  kinetic  energy  equation  will  be  evaluated 
using  Eqs .  (3.2  .9)  to  (3.2.11)  and  only  the  last  terms  in  Eqs.  (  3 .2 . 13)  and 
(3.2. 14)  will  be  calculated . 
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3 .3  Low-Level  Vortlcity  Production 

If  the  complete  vorticity  equation 

(3.3.1)  1  to  +  J )  -- +  dv 

At  3)> 

is  applied  to  a  vorticity  maximum  (  i .  e .,  V(q  +  £)  =  O  ) 
at  a  level  where  CO  is  small,  one  obtains 


(3.3.2)  1  (cj  +  f)  =  -(^4f)V-V+ 

it 

Here  S/it  is  the  local  change  in  a  coordinate  system  moving 
with  the  vorticity  maximum . 

Let  us  now  consider  a  level  within  the  friction  layer  but 
sufficiently  far  from  the  ground  for  the  following  conditions  to  hold: 

0.3.3)  qas  Sj 

and 


(3.3.4)  -  V- V  =  4. 

dp  dp 

Here,  and  CJj,  are  solutions  to  Eqs.  (  3.1.4)  and  (3.1.5). 

Using  Eqs.  (  3.3.3)  and  (  3.3.4),  Eq.  (  3. 3. 2\ may  be 
divided  into  two  components: 
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r  * 

(3-3-5)  L  it 

<Vf)  =  +  *,7xF 

(3.3.6) 

[&<*•')]„  ■ 

where  the  subscripts  Cl 

processes,  respectively 

in  Section  4.6. 

and  h  refer  to  dry-  and  saturated-  adiabatic 
.  The  terms  in  the  above  equations  will  be  computed 

A  method  of  estimating  the  order  of  magnitude  of  the  frictional 

term  in  Eq .  (  3 .3 .5)  will  now  be  described .  The  frictional  force  will  be  assumed 

to  be  expressible  in  the  form 

(3.3.7) 

F  =  -bV 

where  Vb  =  0 

.  We  then  have 

!?•  V*F  =  —  b 

and,  using  Eq.  (  3.3.3), 

• 

(3.3.8) 

l-VxF  =  -b<Jj 
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To  determine  IF  ,  consider  a  column  of  unit  area  extending  from  the  ground 
to  a  height  of  1  km .  The  magnitude  of  the  stress  on  the  bottom  surface  of  the 
column  is  assumed  to  be  given  by 

(3.3.9)  F0=,S,C\/02 

and  the  stress  on  the  top  surface  will  be  neglected .  The  average  force  per 
unit  mass  is  then  approximately 


(3.3.10) 


_?  -W  1 

F  =  (»o  c™  )  cv0 


which  will  be  assumed  to  be  also  the  magnitude  of  IF  at  the  level  in 
question.  Comparing  Eqs.  (  3.3.7)  and  (  3.3.10),  it  will  be  seen  that 


(3.3.11) 


b  = 


(|0  )  C  V, 


z 


V 


Combining  Eqs .  (3.3.8)  and  (3.3.11),  we  then  obtain 


(3.3.12)  Ifc.VxIT  =  -  ('°  ')CV,<,\ 

v  5 

The  fact  that  a  linear  function  of  the  velocity  has  been  used  in  Eq .  (3.3.7)  while 
a  quadratic  one  is  used  in  Eq .  (  3  .3 . 10)  is  not  as  inconsistent  as  would  appear  on 
first  sight.  We  have  only  required  that  Vb  =  °  in  Eq .  (  3 .3 .7)  and 

if  W  r  O  also,  we  can  set 


b  =  b'  v 
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Eq .  (3.3 .7)  then  becomes 


F  =  -  k'  W 


It  is  interesting  to  note  that  it  is  possible  to  determine  the  effect 
of  released  latent  heat  on  Petterssen's  development  equation  (  Petterssen,  1955). 
Integrating  Eq.  (  3.1.2)  from  Jo  ,  the  level  of  non-divergence,  to  p0  we 
have 


i  ft  M  & 


f  Cv5*7s  5f  -  L * 


(3.3.13) 


-  r  a  h 

jp  Cpf 


where  Z*  and  2^,  are  the  heights  of  the  pressure  surfaces  p0  and 

••  <2 
p  respectively.  Applying  the  operator  (l/f)v  to  Eq .  (3.3 . 13)  and 

setting 


ft  |v?p  =-h-v(vf)] 


we  obtain 


1  f  vlZp.  =  -  [vrv<Vf>]p -  ?  r*  * 


(3.3.14) 
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which  is  Petterssen's  development  equation  with  the  geostrophic 
approximation  used  to  represent  the  vorticity  and  the  advective  wind. 


It  is  customary  to  treat  H  in  Eq .  (  3 .3 . 14)  as  heat 
introduced  from  external  sources  and,  in  cases  of  saturated  ascent, 
deal  with  released  latent  heat  by  modifying  y2  e  u>  in  Eq.  (  3.3.14) 
to 


f  ,  R  Lv  /dr  \  1 

V  ff  ■  - —  -  jr  U> 

l  Crf  \+rs  Uf>  4 


Here,  the  subscript  S  denotes  the  saturated  state .  The  term 
(R/C^p)  0-v/(l+rs))  (dr/dp)s  is  equivalent  to  the  right  side 
of  Eq .  (3.1.9).  However,  external  heat  sources  will  be  neglected  here 
and  H  will  be  interpreted  as  released  heat  of  vaporization. 


Writing  OJ  and  making  use  of  Eq.  (3.1.5), 

Eq .  (  3 .3 . 14)  may  be  written 


i.  i 

St  f 


(3.3.15) 


The  last  term  in  Eq .  (3.3. 15)  represents  the  effect  of  released  latent  heat 
and,  had  been  used  on  the  right  side  of  Eq  .  (  3 . 1 . 1),  this 

term  would  be  +  -P )(  . 
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Substituting  Eq .  (  3 . 1 .4)  in  the  third  term  on  the  right  side  of 
Eq.  (  3.3.15),  one  obtains 


(3.3.16) 


If  the  vorticity  advection  at  1000  mb  is  small,  it  is  evident  that  the  effects  of 
dry-  and  saturated -adiabatic  processes  are  proportional  to  ( /dp)  po 
and  (dt^/dp)  po  ,  respectively.  This  result  is  similar  to  that  expressed 
by  Eqs.  (  3.3.5)  and  (3.3.6). 
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Substituting  Eq.  (  3 . 1 .4)  in  the  third  term  on  the  right  side  of 
Eq .  (  3 .3 . 15),  one  obtains 


(3.3.16) 


rt  7  --  -  *4' ^  W 


i>p/p* 


If  the  vorticity  advection  at  1000  mb  is  small,  it  is  evident  that  the  effects  of 
dry-  and  saturated -adiabatic  processes  are  proportional  to  (  p 

and  (du)K/af>)po  ,  respectively.  This  result  is  similar  to  that  expressed 
by  Eqs .  (3.3.5)  and  (3.3.6). 


4. 


A  CASE  STUDY 


4 . 1  General  Synopsis 


The  case  selected  for  study  was  a  cyclone  which  developed  rapidly 
over  the  central  United  States  on  January  21-22,  1959  (Fig.  4.1.1).  The 
system  was  associated  with  an  intense  baroclinic  zone  and  was  accompanied  by 
widespread  heavy  precipitation. 

During  the  period  0000-1200  GMT  of  the  21st  the  frontal  low  moved 
from  western  Texas  to  southeastern  Missouri,  becoming  more  organized  but 
filling  slightly .  Shortly  after  1200  GMT  the  system  began  to  deepen  at  the  rate 
of  approximately  1  mb  hr'  ^ .  By  0000  GMT  of  the  22nc*  it  had  reached  southern 
Ontario  and  had  then  developed  into  a  major  cyclone.  From  the  500  mb  map 
in  the  lower  right  of  Fig .  4 . 1 .1,  it  will  be  seen  that,  in  accordance  with 
Petterssen's  development  hypothesis  (  Petterssen,  1955),  sudden  development 
began  when  the  area  of  appreciable  vorticity  advection  in  mid-troposphere 
spread  over  the  low-level  frontal  zone. 

An  investigation  of  the  dynamic  and  thermodynamic  processes 
occurring  in  the  atmosphere  at  the  onset  of  sudden  cyclogeneis  should  provide  a 
better  understanding  of  the  complex  phenomenon  of  cyclone  development.  It 
was  therefore  decided  to  perform  the  computations  proposed  in  Section  3  using 
the  data  for  1200  GMT  of  the  21st . 

Precipitation  rates  at  1200  GMT  (  Fig.  4.1.2)  were  obtained  from  the 
observed  four -hour  amounts  ending  at  1400  GMT  which  were  then  averaged  over 
areas  centered  at  grid  points  (  see  Fig .  4.3.1)  and  converted  to  hourly 
rates  .  The  time  and  area  averaging  procedure  resulted  in  the  smooth  pattern 
shown  in  Fig.  4.1.2.  Note  the  elongation  of  the  precipitation  pattern  along  the 
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Flg.  4.1.2.  Observed  rate 
of  precipitation .  1200  GMT  Jan . 

21,  1959.  Units:  mm  hr"l. 


frontal  zone .  Furthermore,  an  examination  of  the  surface  data  in  Fig .  4.1.1 
shows  that  most  of  the  precipitation  was  in  the  form  of  rain,  despite  the  fact 
that  this  was  a  mid -winter  situation.  A  possible  explanation  of  this  is  offered 
in  Section  4.4. 

4.2  Kinematically  Computed  Vertical  Velocity 

The  vertical  velocity  was  first  computed  kinematically  at  900,  800, 
600,  400,  300  and  200  mb  by  determining  the  divergence  of  the  mean  wind  between 
successive  levels  ,  For  the  divergence  calculations  a  nomogram  similar  to  the 
one  suggested  by  Fujita  (1955)  was  used . 

Orographic  effects  were  included  in  the  following  manner.  Let 
(  *  »  V  l  be  the  surface  pressure  and  let  p,  be  a  pressure 

level  at  which  k)  is  to  be  determined .  If  dps  /  it  ~  O  , 

OJ(p,)  -  V-V 

a  rps 
V  =  V  $P 
\  1 


where 
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The  above  method  has  been  used  by  a  number  of  Investigators  (  e .  g . ,  Rex,  1958, 
and  Palm^n  and  Holopainen,  1962). 

In  our  computations  Ps  was  put  equal  to  the  standard  atmosphere 
pressure  of  the  area-smoothed  terrain  provided  by  McClain  (1960).  The  pressure 
at  this  smoothed  topography  varies  from  about  1000  to  700  mb  over  the  United 
States,  with  a  large  area  where  ps  >  900  mb .  If  one  used  a  constant  p, 

over  the  entire  region,  the  lowest  convenient  standard  pressure  level  to  be 
selected  for  p,  would  be  600  mb.  However,  if  this  value  of  p,  were 
chosen,  the  vertical  velocity  at  lower  levels  could  not  be  obtained .  It  was  therefore 
decided  to  vary  p,  over  the  region  according  to  the  following  criteria:  • 


ps 

> 

900  mb 

•^3 

II 

900  mb 

900  >  ps 

> 

800  mb 

■y 

II 

800  mb 

Ps 

< 

800  mb 

II 

sL 

600  mb 

The  vertical  velocity  at  higher  levels  was  then  computed  in  the  customary  manner  by 
adding  the  divergence  of  the  mean  wind  of  successive  layers . 

The  fields  of  vertical  motion  at  800,  600,  400  and  200  mb  were  then 
averaged  over  areas  centered  at  grid  points  (  see  Fig.  4.3.1).  The  results  for 
the  region  to  the  east  of  the  Continental  Divide  are  shown  in  Fig.  4.2.1. 

The  accuracy  of  any  kinematic  computation  of  the  vertical  velocity 
normally  decreases  with  height.  Furthermore,  because  of  the  sparsity  of  wind 
data  in  the  region  of  maximum  upward  motion  (see  Fig.  4.1.1),  our  confidence  in 
the  kinematically  computed  values  at  400  mb  and  especially  at  200  mb  is  rather  low. 
Nevertheless,  the  horizontal  projections  of  the  centers  of  maximum  upward  motion 
at  800,  600  and  400  mb  agree  fairly  well  with  the  center  of  heaviest  precipitation 
(see  Fig.  4.1.2).  However,  it  became  evident  when  computing  the  rate  of  release 


Fig .  4.2.1.  Kinematically  computed  vertical  velocity 
1200  GMT  Jan  21,  1959.  Units:  10"S  mb  sec"1. 
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of  latent  heat  (described  In  Section  3.1)  that  In  the  region  of  maximum  pre¬ 
cipitation  the  kinematically  computed  vertical  velocity  accounted  for  only  about 
70  per  cent  of  the  observed  precipitation.  A  possible  reason  for  this  discrepancy 
is  that  the  strong  low-level  convergence  in  this  area  would  indicate  an  appreciable 
import  of  liquid  water  in  the  form  of  stratus  and  stratocumulus  clouds.  According 
to  Cunningham  (1952)  these  clouds  can  contribute  significantly  to  the  precipitation. 

It  is  of  interest  to  note  that  precipitation  occurred  in  Kansas  and 
Oklahoma  even  though  Fig .  4.2.1  shows  downward  motions  at  all  levels  in  this 
region.  This  apparent  inconsistency  is  likely  caused  by  small-scale  turbulence 
or  orographic  effects  which  are  undetectable  in  the  winds  or  smoothed  topography 
used  in  the  computations . 


4.3  Vertical  Velocity  Due  to  Dry-Adiabatic  Processes 

Equation  (3.1.4)  was  solved  for  the  region  shown  in  Fig .  4.3.1 
using  a  vertical  mesh  size  of  200  mb  and  keeping  u)^  =  O  at  p  =  O 
and  jo  =  1000  mb.  The  lateral  boundaries,  on  which  dOa  was  set  equal 

to  zero,  were  located  three  grid  distances  outside  the  boundary  of  the  region 
shown  in  F  ig .  4.3.1. 


In  the  terms  on  the  right  side  of  Eq  .  (3.1 .4),  the  height  data  from 
nine  levels  (  every  100  mb  from  100  to  900  mb  )  were  used .  Let  subscripts  1,  2, 

...,  9  denote  the  levels  p  s  100,  200,  ...,  900  mb.  Then,  in  finite  difference  form, 


/  ^  (<^a)r,+  2  +_  (k-Orv-l  “2. 

'  L  (A  p)  1 


(y..v*?)  g  (V^  .s? 

°  ^  '•n  Ab 
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Kv(vfiL-  M(vf)L 

ap 


where  Ap  =  200  mb,  and  n  takes  on  the  values  2,  4,  6  and  8.  The 

customary  finite  difference  approximations  were  used  for  the  horizontal 
derivatives . 


The  values  of  <T  were  computed  directly  from  the  temperatures 
also  using  the  data  at  nine  levels  from  100  to  900  mb.  In  finite  difference  form. 


(6),  *P 


Equation  (3.1.4)  was  first  solved  in  the  form 


where  0o(  p)  is  the  area -averaged  value  of  <J(Xfy9p) 

The  solution  at  600  mb  is  reproduced  in  Fig.  4.3.2.  Comparing  Figs.  4.3.2 
and  4,2.1,  it  is  evident  that  the  solution  of  Eq  .  (  4 .3 . 1)  grossly  underestimates 
the  magnitudes  of  the  maximum  downward  and  upward  motion,  especially  that  of 
the  latter . 


One  reason  for  this  discrepancy  is  due  to  the  neglect  of  horizontal 
variations  of  C  .  Fig .  4 .3 .3  shows  the  variations  of  /  00  at 
600  mb.  This  quantity  varies  from  less  than  0.5  to  greater  than  2.0.  It  will 
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Fig.  4.3.1.  Region  of 
solution  and  grid  points  used 
for  the  numerical  computations . 


Fig.  4.3.2.  Vertical 
velocity  at  600  mb  associated 
with  dry-adiabatic  processes 
keeping  the  static  stability 
constant.  1200  GMT  Jan.  21, 
1959.  Units:  10'^  mb  sec"^. 


Fig.  4.3,3.  Static  stability 
at  600  mb  divided  by  its  areal 
average  value.  1200  GMT  Jan.  21, 
1959.  * 
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be  noted  also  that,  as  one  might  expect,  the  air  in  the  frontal  zone  is  more 
statically  stable  than  either  the  warm  or  cold  air. 

Equation  (3 . 1 .4)  was  then  solved,  permitting  to  vary  along  the 
isoharic  surfaces  (  Fig .  4 .3 .4) .  Comparing  Figs.  4.3.4  and  4.3.2,  it  will  be 
seen  that  one  effect  of  letting  0  vary  is  to  shift  the  maxima  of  upward  and 
downward  motions  towards  regions  of  lower  stability  and  to  Increase  their 
magnitudes .  This  amplification  is  due  to  the  circumstance  that  O / C0  <  I 
in  these  regions .  Comparing  Figs.  4.3.4  and  4.2.1,  note  that  the  magnitude 
of  the  maximum  downward  motion  is  correctly  computed  from  Eq.  (3.1.4) 
although  the  center  is  too  far  west .  It  is  likely  that  inclusion  of  terrain  effects 
in  the  numerical  solution  would  have  resulted  in  better  agreement .  The  effect 
of  neglecting  orographic  influences  is  especially  noticable  in  the  northwestern 
section  of  the  region .  However,  the  computed  maximum  of  upward  motion,  though 
larger  than  in  the  case  of  constant  static  stability,  is  still  seriously  underestimated . 

4 .4  Effect  of  Release  of  Latent  Heat  on  the  Vertical  Velocity 

To  solve  Eq .  (3.1.5),  the  rate  of  release  of  latent  heat  (  H  ) 
must  be  known .  This  quantity  was  computed  using  the  method  described  in 
Section  3.1;  the  results  are  shown  in  Fig.  4.4.1. 

Equation  (3 . 1 ,5)  was  then  solved,  permitting  to  vary  and  using 
the  same  region  and  boundary  conditions  as  described  in  Section  4.3 .  The 
solution  at  600  mb  is  presented  in  Fig .  4  .4 . 2 .  The  solutions  at  800  and  400  mb 
were  similar  with  maxima  of  about  -3.5  and  -4.0  x  10"^  mb  sec'*,  respectively. 
Even  though  there  was  no  release  of  latent  heat  at  200  mb,  at  this  level 

was  not  zero.  However,  its  magnitude  was  everywhere  less  than  0.3  x  10 "3  mb 
sec'*. 
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Fig.  4.3.4.  Vertical  velocity  associated  with  dry-adiabatic 
processes  permitting  the  static  stability  to  vary .  1200  GMT  Jan . 
21,  1959.  Units:  10‘3  mb  sec"*. 
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Pig.  4.4.1.  Rate  of  release  of  latent  heat  computed  from  the  observed 
precipitation.  1200  GMT  Jan  21,  1959.  Units:  10^  ergs  gm-1  sec-1. 


Fig.  4 .4 .2,  Vertical  velocity  at  600  mb 
associated  with  saturated -adiabatic  processes 
permitting  the  static  stability  to  vary .  1200 
GMT  Jan.  21,  1959.  Units:  10'3  mb  sec-1. 
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From  Fig.  4.4.2  it  will  be  seen  that  the  effect  of  released  latent 
heat  is  to  amplify  the  upward  motion  in  the  region  of  heaviest  precipitation  and 
contribute  to  a  slight  downward  motion  in  the  surrounding  region.  This  same 
effect  has  been  observed  by  Smagorinsky  (1956),  Aubert  (1957)  and  Pedersen  (1962). 
Furthermore,  this  result  is  in  agreement  with  the  theoretical  solution  Eliassen 
(1959)  obtained  in  considering  the  effect  of  an  elliptically-shaped  heat  source  in 
a  baroclinic  zonal  current . 

A  possible  explanation  why  the  precipitation  in  the  case  studied  was 
predominantly  rain  is  the  following.  If  there  exists  initially  a  large  region  of 
relatively  uniform  upward  motion  such  as  that  given  by  Fig.  4.3.4  in  the  region 
of  heaviest  precipitation,  the  vertical  velocity  will  be  amplified  in  those  areas 
which  first  attain  saturation  and  damped  in  the  surrounding  areas .  In  the  case 
under  study,  condensation  presumably  first  occurred  in  the  maritime  tropical 
air  near  the  front,  producing  rain.  Further  north,  where  the  temperatures 
were  low  enough  to  yield  snow,  the  upward  motion  was  decreased  and  the 
attainment  of  saturation  delayed . 

The  solutions  of  Eqs .  (  3 .1 .4)  and  (3 . 1 .5),  permitting  0  to  vary  in 
both  equations,  were  added;  the  results  are  shown  in  Fig.  4.4.3.  Comparing 
Figs.  4.4.3  and  4.2.1  we  see  that  the  agreement  is  satisfactory  at  800  and  600 
mb,  except  that  the  numerical  solutions  still  give  the  maximum  downward 
motion  too  far  west  and  do  not  detect  the  orographically  induced  upward 
motion  in  the  northwest  section  of  the  region.  The  maximum  upward  motion 
is  slightly  underestimated  but  this  would  be  expected  since  frictional  convergence 
has  been  neglected .  The  agreement  is  poorer  at  400  and  more  so  at  200  mb. 
However,  as  has  been  pointed  out  in  Section  4.2,  our  confidence  in  the 
kinematically  computed  vertical  velocities  at  these  levels  is  rather  low . 
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Fig.  4.4.3.  Sum  of  the  vertical  velocities  associated  with 
dry-  and  saturated -adiabatic  processes  permitting  the  static 
stability  to  vary.  1200  GMT  Jan.  21,  1959.  Units:  lO'^mbsec"*. 
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To  sum  up,  it  can  be  stated  that,  in  regions  of  heavy  precipitation, 
one  cannot  compute  the  vertical  velocity  with  adequate  accuracy  by  numerical 
means  unless  one  includes  the  effect  of  released  heat  of  vaporization. 


4.5  Influence  of  Released  Latent  Heat  on  the  Changes  in  Kinetic  Energy 

The  induced  ageostrophic  wind  components  and  were 

determined  from  Eqs.  (  3.2.4)  and  (3.2.5)  at  900,  700,  500  and  300  mb  using  the 
values  of  COa  and  obtained  from  the  solutions  of  Eqs.  (  3.1.4)  and 

(3.1.5).  The  results  at  900  and  300  mb  are  shown  in  Fig .  4.5.1.  Note  the 
reversal  in  direction  of  both  and  between  the  two  levels . 

The  right  hand  sides  of  Eqs.  (  3.2.9)  to  (  3.2.11)  were  then 
evaluated  over  the  region  given  in  Fig.  4.3.1  for  four  layers,  each  200  mb 
deep,  between  1000  and  200  mb.  Table  4.5.1  shows  the  results  of  computing 
these  quantities  and  adding  them  to  obtain 


The  table  also  shows  the  amount  of  kinetic  energy,  the  time  required  for  the 
kinetic  energy  to  double  (  assuming  K  4  /  it  to  be  constant),  and 
the  total  rate  of  release  of  latent  heat  within  each  layer  (  H  ^  ).  For 
convenience,  the  layers  are  identified  by  numbers,  thus:  1  (1000  to  800  mb); 
2  (  800  to  600  mb);  3  (  600  to  400  mb);  and  4^400  to  200  mb) . 
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Table  4.5.1.  Comparison  of  the  terms  in  Bq .  (4.5.1). 
For  definition  of  symbols,  see  text. 


Term 

Units 

Layer 

Totals 

1 

2 

3 

4 

Ht 

.„25 

10  ergs 

1.95 

3.50 

9.57 

18.44 

33.46 

,„20 

10  ergs 

-1 

sec 

27.9 

77.1 

86.5 

37.2 

228.7 

[aXj/w], 

,„20 

10  ergs 

-1 

sec 

0.34 

0.57 

1.23 

-0.01 

2.13 

,_20 

10  ergs 

-1 

sec 

1.26 

1.09 

3.01 

3.89 

9.25 

0Kf/Jt]h 

ir.20 

10  ergs 

-1 

sec 

1.28 

-0.39 

0.54 

3.05 

4.48 

d  Ky/dt 

lri20 

10  ergs 

-1 

sec 

2.88 

1.27 

4.78 

6.93 

15.86 

Doubling 

time 

hours 

18.8 

76.4 

55.7 

73.7 

58.4 

The  effect  of  neglecting  surface  friction  is  reflected  in  the  small 
doubling  times  computed  for  layer  1.  However,  the  doubling  times  computed  for 
layers  2,  3  and  4  appear  reasonable,  although  there  is  some  variation  which  can 
be  attributed  to  one  or  more  of  the  following  causes: 

(i)  horizontal  import  of  kinetic  energy  by  the  ageostrophic  wind; 

(ii)  rate  of  working  by  the  pressure  due  to  the  rotational  component 
of  the  ageostrophic  wind;  and 

(iii)  errors  arising  from  inaccuracy  in  the  input  data . 


Fig.  4.5.1.  Irrotatlonal  ageostrophic  wind  components  at  900  and  300  mb. 
The  charts  on  the  left  and  right  show  the  components  associated  with  dry-  and 
saturated-adiabatic  processes,  respectively .  Solid  lines:  streamlines .  Dashed 
lines:  Isotachs  In  units  of  m  sec'*.  1200  GNffT  Jan.  21,  1959. 
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It  may  be  noted  also  that  in  each  layer  H  ^  is  an  order  of  magnitude 
larger  than  5  K  4  /  3 1  .  Comparing  /dt]  and 

C  »tq  /»t  ]  h  we  see  that  in  layers  1  and  4,  both  terms  are  approximately 

equal,  although  at  mid-levels  [»K$  /  at]  K  is  rather  small  and  is 

actually  negative  in  layer  .  . 


Further  insight  is  gained  if  the  individual  terms  in  Eqs.  (3.2.10)  and 
(3.2.11)  are  examined.  Let  us  write  these  equations  symbolically  as 


(4.5.2) 


(4.5.3) 


[  |p]4  -  I„  .  W. 

[if'],  ■  u  * 


where  Tyo.  and  represent  the  net  vertical  import  of  kinetic  energy 

by  6Ja  and  respectively,  and  Wfl  and  denote  the  rate  of 

working  by  the  pressure  by  and  .  The  results  are  shown  in  Table 

4.5.2. 


Table  4.5.2.  Comparison  of  the  terms  in  Eqs.  (  4.5.2)  and  (4.5.3). 
For  definition  of  symbols, see  text.  Units:  1020  ergs  sec'*. 


Term 

Layer 

Totals 

1 

2 

3 

4 

Iva 

-0.59 

-0.66 

-0.03 

0.43 

-0.86 

Wa 

1.85 

1.75 

3.04 

3.47 

10.11 

IvK 

-0.47 

-0.88 

-0.18 

0.87 

-0.29 

Wh 

1.75 

0.49 

0.36 

2.18 

4.78 
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The  terms  Xya  and  XVK  are  negative  In  the  lower  troposphere 
and  positive  at  higher  levels,  indicating  a  net  transport  of  kinetic  energy  from  the 
lower  to  upper  troposphere .  Moreover,  they  are  generally  of  smaller  magnitude  than 
the  terms  and  VVj,  •  These  latter  terms  are  positive  at  all  levels.  The 

mid -troposphere  minimum  of  [aK^/at]  h  ,  noted  above,  evidently  results 
from  Wh  .  The  vertical  variation  of  this  last  term  can  be  readily  explained. 


At  mid-levels,  where  the  rate  of  release  of  latent  heat  is  a  maximum, 
/dp  is  near  zero.  Consequently  the  induced  ageostrophic  velocity 
is  small,  and  therefore  so  is  W  ^  • 


At  low  levels,  the  center  of  maximum  convergence  caused  by  released 
latent  heat  is  located  in  a  region  of  relatively  low  contour  heights .  The  induced  inflow 
(  see  Fig.  4 .5 . 1)  towards  this  center  is  therefore  generally  in  the  direction  of 
decreasing  pressure,  resulting  in  positive  work  being  done  on  the  air  parcels . 


At  high  levels,  the  center  of  maximum  divergence  is  situated  in  a  region 
of  relatively  high  contour  heights.  This  is  because  in  the  upper  troposphere  a  warm 
ridge  lies  just  east  of  the  location  of  the  surface  frontal  low .  Hence  the  induced  out¬ 
flow  is  again  generally  in  the  direction  of  decreasing  pressure. 

The  term  K  has  a  similar,  though  less  pronounced,  vertical 

variation . 


Although  only  one  case  has  been  investigated  here,  it  is  likely  that 
similar  results  would  be  obtained  from  further  studies .  It  is  typical  of  frontal  systems 
that  precipitation  occurs  in  a  region  of  relatively  low  contour  heights  in  the  lower 
troposphere  surmounted  by  a  warm  ridge  at  higher  levels .  Therefore  it  is  to  be 
expected  that  in  such  systems  accompanied  by  moderate  or  heavy  precipitation, 
released  latent  heat  makes  a  significant  contribution  to  the  kinetic  energy  changes  at 
low  and  high  levels  but  has  little  effect  in  mid-troposphere . 
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Next,  the  energy  conversion  integrals 


(4.5.4)  <V 
and 

(4.5.5)  Ch--JJ“h«iAy 


were  computed .  The  results  are  shown  in  Table  4.5.3. 


Table  4.5.3.  Comparison  of  Eqs .  (  4  .5 .4)  and 
(4.5.5).  Units:  l(r®  ergs  sec”*. 


Term 

1 - 

Layer 

Totals 

1 

2 

3 

4 

c* 

0.46 

8.28 

25.92 

37 

.63 

72.29 

Ch 

8.76 

27.71 

40.77 

33 

.15 

110.39 

Except  at  level  4,  the  contribution  of  released  latent  heat  to  the 
energy  conversion  processes  exceeds  that  caused  by  dry-adiabatic  motion. 
Comparing  the  values  in  Table  4.5.3  with  the  computed  magnitudes  of 

and  given  in  Table  4 .5 . 1,  it  is 

evident  that  the  other  terms  in  Eqs .  (  3.2.13)  and  (  3 .2 . 14)  (not  computed)  must 
in  general  be  of  opposite  sign  to  Ca  and  C), 
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4.6  The  Role  of  Released  Latent  Heat  in  Vortlcity  Production 

The  change  in  absolute  vortlcity  following  the  motion  of  the  sea- 

level  center  (  see  Fig.  4.1.1)  between  0000  GMT  of  the  21st  and  0000  GMT  of 
nd 

the  22  was  observed  to  be 

(4.6.1)  Lfa.  +  t)  -  0.88 X  10"*  sec'1  <24  hr)'1 

it  ’ 


Since  the  rate  of  change  of  vorticlty  was  fairly  uniform  throughout  the  24-hour 
period,  the  value  given  by  Eq .  (4.6.1)  may  be  taken  as  the  instantaneous  rate 
of  change  at  1200  GMT  of  the  21st. 


The  right  hand  sides  of  Eqs  .(3.3 .5)  and  (3.3 .6)  were  computed 

st 

from  conditions  at  1200  GMT  of  the  21  ,  using  the  following  numerical  values: 


V 

-f  - 
c  = 
\4  = 


1.79xl0'4  sec'1 

V  = 

20  m  sec  1 

0.87  x  10  4  sec  1 

huia/dp  — 

0.4  x  10  ^  sec  1 

0.005 

d^/dp  = 

1.5  x  10  **  sec  1 

-1 

10  m  sec 

The  numerically  computed  convergence 


ScOfl,  —  1.9x10^  sec  1 

^  *P 

agreed  well  with  the  value  determined  directly  from  the  wind  observations 
(2.3  x  10  5  sec  V  The  following  results  were  obtained: 

W  'Vx  F  =  -3.87  x  lo”4  sec"1  (24  hr) 
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sec  1  (24  hr)  * 
sec  1  (24  hr)  1 

sec  1  (24  hr)  1 

In  view  of  the  approximations  made  in  computing  the  term 
ft?  *  V  x  IF  ,  the  agreement  between  the  values  given  by  Eqs.  (4.6.2) 

and  (  4 .6 . 1)  must  be  considered  fortuitous .  Nevertheless,  the  frictional  term 
must  be  large  to  compensate  for  the  tendency  for  rapid  intensification  resulting 
from  convergence.  Thus  the  observed  changes  in  vorticity  evidently  resulted 
from  a  small  imbalance  between  two  large  quantities,  production  of  vorticity 
by  convergence  and  destruction  by  friction . 

If  condensation  suddenly  ceased  in  the  system, 

=  -2.95  x  10"4  sec'1  (24  hr)"1 

which  would  correspond  to  rapid  weakening.  Conceivably  this  could  be 
one  reason,  along  with  an  increased  drag  coefficient,  for  the  rapid  decay 
observed  when  hurricanes  enter  land  and  are  cut  off  from  their  moisture  supply. 
It  is  not  intended  to  imply  that  the  quasi-geostrophic  equations  (3.1,4)  and 
(3.1.5)  are  necessarily  adequate  in  the  case  of  a  hurricane  to  depict  the  fields 


(e^  =  0.92x10 

+  f)(<Sulh/af>)  =  3.45x10 


-4 

4 


Thus 


(4.6.2)  [ftty+V]  =  °-50x  10"4 
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of  vertical  motion  and  convergence .  Nevertheless,  in  a  hurricane  one  would 
expect  released  latent  heat  to  have  an  effect  on  these  fields  similar  to  that 
demonstrated  in  this  paper . 

Since  even  small  values  of  low-level  convergence  (of  the  order  of 
0.5  to  1.0  x  10  sec  *)  would  by  themselves  produce  vorticity  changes  larger 
than  those  observed  in  the  atmosphere,  it  is  likely  that  at  all  stages  in  the 
evolution  of  a  low  the  opposing  effects  of  friction  and  convergence  very  nearly 
balance.  The  initiation  of  condensation  would  then  result  in  an  imbalance  and  an 
intensification  of  the  system  with  an  accompanying  change  in  the  frictional  effect, 
caused  by  a  modification  of  the  velocity  field,  restoring  the  state  of  near-balance . 


5. 


CONCLUDING  REMARKS 


Although  only  one  synoptic  situation  has  been  examined  here,  the 
results  appear  general  enough  to  permit  the  tentative  conclusion  that  they  are 
typical  of  lows  associated  with  strong  baroclinity  and  accompanied  by  moderate 
to  heavy  precipitation .  The  results  for  the  case  investigated  may  be  summarized 
as  follows . 


In  the  case  studied,  the  conventional  method  of  obtaining  the 
vertical  velocity  from  the  diagnositc  OJ  -equation  by  assuming  constant  static 
stability  and  dry-adiabatic  motion  proved  to  be  inadequate .  The  maximum 
upward  motion  obtained  in  this  way  was  only  about  one -fourth  of  the  value 
computed  by  the  kinematic  method .  Oily  by  permitting  the  static  stability  to 
vary  and  by  including  the  influence  of  released  latent  heat  ( the  latter  modification 
being  the  more  important)  was  satisfactory  agreement  between  the  numerically 
and  kinematically  computed  values  obtained. 

Accompanying  the  amplifying  effect  that  released  heat  of 
vaporization  has  on  the  vertical  velocity  in  mid-troposphere  was  an  increase 
in  low-level  convergence  and  high-level  divergence.  Furthermore,  owing  to 
the  circumstance  that  the  precipitation  was  associated  with  a  region  of 
relatively  low  contour  heights  in  the  lower  troposphere  surmounted  by  a  warm 
ridge  at  higher  levels,  the  induced  ageostrophic  horizontal  motion  was  in  the 
direction  of  decreasing  pressure  at  both  low  and  high  levels,  resulting  in 
production  of  kinetic  energy.  This  production  was  of  the  same  order  of 
magnitude  as  that  caused  by  dry-adiabatic  processes .  At  middle  levels, 
released  latent  heat  had  little  effect  on  the  changes  in  kinetic  energy . 
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The  maximum  low-level  convergence  associated  with  release  of  latent 
heat  was  about  three  times  as  large  as  that  arising  from  dry-adiabatic  processes, 
excluding  frictional  convergence .  One  might  therefore  conclude  that  the  pro¬ 
duction  of  vorticity  was  closely  related  to  the  rate  of  release  of  latent  heat.  Since 
the  vorticity  production  predicted  by  the  numerically  computed  convergence,  which 
agreed  well  with  the  value  calculated  directly  from  the  observed  winds,  greatly 
exceeded  the  observed  production,  one  might  further  conclude  that  this  production 
was  very  nearly  balanced  by  frictional  effects . 

It  should  be  pointed  out  that  the  modifying  effects  of  released  heat 
of  vaporization  can  come  into  play  only  after  development  has  begun  and  con¬ 
densation  has  been  achieved .  Therefore  one  cannot  say  that  released  latent 
heat  initiates  cyclones  but  rather  that  it  affects  their  subsequent  growth . 
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